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We show that electronic-entropy effects in the size- 
evolutionary patterns ol relatively small (as small as 20 
atoms), simple-metal clusters become prominent already at 
moderate temperatures. Detailed agreement between our 
finite-temperature-shell-correction-method calculations and 
experimental results is obtained for certain temperatures. 
This agreement includes a size-dependent smearing out of 
fine-structure features, accompanied by a measurable reduc- 
tion of the heights of the steps marking major-shell and sub- 
shell closings, thus allowing for a quantitative analysis of clus- 
ter temperatures. 



PACS numbers: 36.40.-c, 65.90.-H, 05.30.Fk 

Since the discovery |jj of electronic-shell-structure fea- 
tures in the abundance spectra of sodium clusters, sim- 
ilar features (the major ones corresponding to the de- 
generacies of a spherically symmetric mean-field po- 
tential \$Q) have been routinely observed M in the 
size-evolutionary patterns (SEP's) associated with other 
single-particle properties of both alkali- and noble-metal 
clusters. Specifically, such patterns pertain to ioniza- 
tion potentials (IP's) |-0], electron affinities (EA's) ||,|, 
monomer separation energies (MSE's) EO], and fission 
dissociation energies | |L1|. 

It was early realized |12(| that the secondary features in 
the mass spectra required consideration of deformed clus- 
ter shapes JqJia] . When triaxial (ellipsoidal) shapes were 
considered in the framework of Shell Correction Methods 
(SCM) [OTL61, substantial overall systematic agreement 
was achieved |lj,p5[ between theory and experimental 
observations pertaining to the major and the fine struc- 
ture of the aforementioned SEP's. 

While deformation effects have been extensively stud- 
ied, an understanding of the physical origins of thermal 
effects and their relation to the SEP's is still lacking, even 
though the experiments are necessarily made with clus- 
ters at finite temperatures. Moreover, experimental de- 
termination of cluster temperatures remains a challenge, 
motivating the development of theoretical methods capa- 
ble of extracting such information. 

While thermodynamic entropic contributions associ- 
ated with the ionic degrees of freedom can be ob- 
tained from first-principles molecular-dynamics simula- 
tions ||17|Jlq 1 , or from considerations of shape fluctuations 
in simpler models |19]-|2l||, for sufficiently large simple- 
metal clusters M^- with N > 20, the electronic entropy 



(which has not as yet been included in such studies) dom- 
inates the thermal characteristics of the shell structure, 
even at moderate temperatures. The prominent thermal 
effects associated with the electronic entropy are the fo- 
cus of this paper. 

Without the simultaneous consideration of shape fluc- 
tuations, thermal effects pertaining to the electronic de- 
grees of freedom and their relevance for the understand- 
ing of mass abundance spectra have been considered in 
the case of spherical neutral sodium clusters J22| , and in 
a recent report 23 on the thermodynamics of neutral 
sodium clusters with axially symmetric shapes. In con- 
trast to our findings, these studies have suggested that 
electronic-entropy effects at moderate temperatures are 
not important for clusters with less than 100 atoms. 

The theoretical method used in this paper is a finite- 
temperature (FT) -SCM developed by us, which incorpo- 
rates all three of the aforementioned aspects, namely, tri- 
axial deformations, entropy of the electrons, and thermal 
effects originating from shape fluctuations. Furthermore, 
through a direct comparison with experimental measure- 
ments, we demonstrate that this method can be employed 
for determining cluster temperatures. 

Since the number of delocalized valence electrons is 
fixed for a given cluster, M^ , we need to use P4,p5[ the 
canonical ensemble in calculating their thermodynamical 
properties. For determining the free energy, F(j3, N, x) 
(f'3 = 1/T), which incorporates the electronic entropy, we 
separate it, in analogy with the zero-temperature limit 
p4| , into a smooth, liquid-drop-model part, -Fldm, and 
a Strutinsky-type |2(| shell-correction term, AF s h. The 
shell correction term is specified as the difference AF s h = 
F sp — F sp , where F sp is the canonical free energy of the 
valence electrons viewed as independent single particles 
in their effective mean-field potential. For calculating 
the canonical F sp , we adopt a number-projection method 



according to which 
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where N e is the number of electrons and fi is the chemi- 
cal potential of the equivalent grand-canonical ensemble. 
The grand-canonical partition function, Zqq, is given by 

Z GC ((3,P^)=l[(l + e-^^), 



where e% are the single-particle levels of the modified Nils- 
son hamiltonian pertaining to triaxial shapes El], The 



temperature-dependent average F sp |£8[ is specified us- 
ing the same expressions as for _F sp , but with a set of 
smooth levels || defined as e % = E™ c (i) - E™ c (i - 1), 

where E°p C (N e ) is the zero-temperature Strutinsky av- 
erage of the single-particle spectrum of an anisotropic, 
triaxial oscillator (see section II. C. of Ref. fl4[| ). 

The LDM term -Fldm consists of three contributions; 
a volume, a surface, and a curvature term. Since vol- 
ume conservation during deformation is assumed, we 
need not consider the temperature dependence of the 
corresponding term when calculating observables, such 
as IP's and MSE's, associated with processes which do 
not change the total number of atoms, N. The exper- 
imental temperature dependence of the surface tension, 
a(T) = cq — c\(T — T m p)) is taken from standard tables 
p9[ (T mp are melting-point temperatures Q), but nor- 
malized to yield the a(T = 0) value used in our earlier 
calculations G3]. Since no experimental information is 
available regarding the curvature coefficient, A c , we as- 
sume the same relative temperature dependence as for 
<j{T) normalized to the T = value used earlier E3. Fi- 
nally, the expansion of the Wigner-Seitz radius due to the 
temperature is determined according to the cofficient of 
linear thermal expansion [go). With these modifications, 
the remaining steps in the calculation of Fldm follow 
closely section II. A. of Ref. jlij (for simplicity, the work 
function W is assumed temperature independent). 

According to the general theory of thermal fluctuating 
phenomena |3l|, the triaxial shapes of the clusters, speci- 
fied by the /3h and 'Jh Hill- Wheeler parameters |32|] , will 
explore Q the free-energy surface, F(0,N,x;Ph,7h), 
with a probability, 

P(Ph,1h) = Z- 1 vv[-PF{P,N,x;Ph,1h)\ , 

the quantity Z being the classical Boltzmann-type par- 
tition function, 



dTexp[-/3F(/3,Ar,:r;/3 ff ,7tf)] 



and dr = (3jj\ sm(3-fH)\df3iid'yH the proper volume ele- 
ment J34J. Thus, finally, the free energy, (F(/3,N,x)), 
averaged over the shape fluctuations can be written as 

(F(/3,N,x)) = J dTF(/3,N,x;l3 H , lH )P(f3 H ,~f H ) . 

We will present results pertaining to IP's and MSE's 
(Mjy — > Mjy_ 1 + M), which at finite temperatures are 
defined as, 

I N = (F(0, N,x = +1)) - (F(/3, N,x = 0)) , 

and 

D U = < F ^< N-l,x = +1)) + (F(/3, N=l,x= 0)> 
-{F(/3,N,x = +l)), 



respectively. 

Our calculations (solid dots) for the IP's of Kjv clusters 
for three temperatures, T — 10 K, 300 K, and 500 K, are 
displayed in Fig. 1, and are compared with the experi- 
mental measurements H (open squares; the experimen- 
tal uncertainties are 0.06 eV for N < 30 and 0.03 eV for 
N > 30). As was the case with our earlier T = K results 
|lj|, the T = 10 K theoretical results exhibit the follow- 
ing two characteristics: (i) Above N — 21, a pronounced 
fine structure between major-shell closures which is not 
present in the experimental measurements; (ii) Steps at 
the major-shell closures which are much larger than the 
experimental ones |3q] (three-to-five times for N = 40, 
58, and 92, and two-to-three times for N = 8 and 20). 

The agreement between theory and experiment is sig- 
nificantly improved at T = 300 K. Indeed, in comparison 
with the lower-temperature calculations, the T = 300 
K results exhibit the following remarkable changes: (i) 
Above N = 21, the previously sharp fine-structure fea- 
tures are smeared out, and as a result, the theoretical 
curve follows closely the smooth modulations of the ex- 
perimental profile. In the size range N = 21 — 34, three 
rounded, hump-like formations (ending to the right at 
the subshell closures at N = 26, 30, and 34) survive in 
very good agreement with the experiment (the sizes of 
the drops at N = 26, 30 and 34 are comparable to the 
experimental ones |39] ) ; (ii) The sizes of the IP drops at 
N = 20, 40, 58, and 92 are reduced drastically and are 
now comparable to the experimental ones. In the size 
range N < 20, the modifications are not as dramatic. 
Indeed, one can clearly see that the pattern of odd-even 
alternations remains well defined, but with a moderate 
attenuation in amplitude, again in excellent agreement 
with the experimental observation. 

For T = 500 K, the smearing out of the shell structure 
progresses even further, obliterating the agreement be- 
tween theory and experiment. Specifically, the steps at 
the subshell closures at N = 26 and 30, as well as at the 
major-shell closures at N = 40, 58, and 92 are rounded 
and smeared out over several clusters (an analogous be- 
havior has been observed in the logarithmic abundance 
spectra of hot, singly cationic, copper, silver, and gold 
clusters |3ffl). At the same time, however, the odd-even 
alternation remains well defined for N < 8. We fur- 
ther notice that, while some residue of fine structure sur- 
vives in the range N — 9 — 15, the odd-even alternations 
there are essentially absent (certain experimental mea- 
surements [B8| of the IP's of hot Na^r clusters appear to 
conform to this trend). 

To ascertain the relative weight of the two thermal pro- 
cesses incorporated in our FT-SCM, we display in the 
middle panel of Fig. 1 (uppermost curve with open cir- 
cles) the theoretical IP's at T = 300 K in the case when 
the electronic entropy is neglected. A comparison with 
the results (solid dots) when both electronic and shape- 
fluctuations entropic contributions are included demon- 



strates the principal role of the electronic entropy in shap- 
ing the thermal effects of the SEP's. 

Fig. 2 displays for two temperatures (10 K and 300 K) 
the FT-SCM results (solid dots) for the MSE's associated 
with the K% clusters, along with available experimental 
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measurements 1 10 (open squares) in the size range N — 
4—23. Compared to the T = 10 K results, the theoretical 
results at T — 300 K are in better agreement with the 
experimental ones due to an attenuation of the amplitude 
of the alternations (e.g., notice the favorable reduction in 
the size of the drops at N = 9, 15, and 21). In spite of this 
amplitude attenuation, it is remarkable that the T = 300 
K SCM results in this size range preserve in detail the 
same relative pattern as the T = 10 K ones (in particular, 
the well-defined odd-even oscillations in the range N = 
4—15 and the ascending quartet at N = 16 — 19 followed 
by a dip at N = 20). 

As a last example, Fig. 3 displays for three temper- 
atures (10 K, 800 K, and 2000 K) our FT-SCM results 
for the IP's of Ag^r clusters, and compares them with 
available experimental measurements [p 39 1. Again, go- 
ing from the T = 10 K to the T = 800 results, we ob- 
serve that an attenuation of the amplitude of alternations 
brings theory and experiment in better agreement (e.g., 
in the latter case, the sizes of the theoretical IP drops at 
TV = 6, 8, 14, 20, 26, 30, and 34 are comparable to the 
experimental ones). Finally, at T = 2000 K, both major 
and fine structure tend to vanish for N > 8. 

In conclusion, we showed that the SEP's of the single- 
particle properties of simple-metal clusters are governed 
by the interplay of ellipsoidal deformations and ther- 
modynamics (entropy) of the electronic degrees of free- 
dom, while the entropic contribution of shape fluctua- 
tions plays a smaller role E0]. We further demonstrated 
that electronic-entropy effects are reflected in prominent 
experimental signatures already at moderate tempera- 
tures and for relatively small sizes. This behavior, which 
had not been previously understood in the case of metal 
clusters, correlates with an order-of-magnitude estimate 
J54| from theoretical nuclear physics. Accordingly, the 
shell-correction term in the free energy depends on T as 
AF sh (T) = AF sh (0)r/sinhT, which for r > 1 takes the 
form 2AF sh (0)rexp(-r) with r = 2-K 2 T/(huj sh ), with 
fiujsh being the energy spacing between shells. Because of 
the large 2ir 2 factor, the shell-structure effects decrease 
rapidly, even for temperatures that are only a fraction 
of the shell spacing (see p. 608 of Ref. [||). Thus, 
the shell structure "melts" for temperatures as low as 
T/ihuish) = 0.25 Q- The FT-SCM calculations pre- 
sented here correspond to temperature ranges well below 
this ratio, where analysis of the extent of attenuation of 
shell-structure signatures allows determination of cluster 
temperatures. 

Compared to nuclear and atomic physics, which have 
heretofore provided the prototypes for fermionic shell 
structure, the SEP's of metal clusters stand apart, since 



the experimentally available SEP's of nuclei (e.g., those 
of neutron separation energies p2y ) and of atoms (i.e., 
IP's p2[ ) correspond strictly to zero-temperature. Fi- 
nally, it is natural to conjecture that thermal (in partic- 
ular, electronic-entropy) effects will influence the height 
of fission barriers and fragmentation channels [E3[ . 
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FIG. 1. IP's of Kjv clusters at three temperatures, T — 10 
K, 300 K, and 500 K. Solid dots: Theoretical FT-SCM re- 
sults. Open squares: experimental measurements H. The 
uppermost curve (open circles) in the middle panel displays 
theoretical results when the electronic entropy is neglected. 
The y-scale in this instance is to the right. 



FIG. 2. MSE's of K^ clusters at two temperatures, 
T = 10 K, and 300 K. Solid dots: Theoretical FT-SCM re- 
sults. Open squares: experimental measurements [|L0|. To 
facilitate comparison, the SCM results at the higher tem- 
perature have been shifted by 0.07 eV, so that the theoret- 
ical curves at both temperatures refer to the same point at 
N = 10. 



FIG. 3. IP's of Agjv clusters at three temperatures, 
T = 10 K, 800 K, and 2000 K. Solid dots: Theoretical 
FT-SCM results. Open squares: experimental measurements 
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